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Abstract 

Here we analyze synaptic transmission from an information-theoretic 
perspective. We derive closed-form expressions for the lower-bounds on 
the capacity of a simple model of a cortical synapse under two explicit 
coding paradigms. Under the "signal estimation" paradigm, we assume 
the signal to be encoded in the mean firing rate of a Poisson neuron. The 
performance of an optimal linear estimator of the signal then provides 
a lower bound on the capacity for signal estimation. Under the "signal 
detection" paradigm, the presence or absence of the signal has to be de- 
tected. Performance of the optimal spike detector allows us to compute 
a lower bound on the capacity for signal detection. We find that single 
synapses (for empirically measured parameter values) transmit informa- 
tion poorly but significant improvement can be achieved with a small 
amount of redundancy. 

1 Introduction 

Tools from estimation and information theory have recently been applied by researchers 
(Bialek et. al, 1991) to quantify how well neurons transmit information about their random 
inputs in their spike outputs. In these approaches, the neuron is treated like a black-box, 
characterized empirically by a set of input-output records. This ignores the specific nature 
of neuronal processing in terms of its known biophysical properties. However, a systematic 
study of processing at various stages in a biophysically faithful model of a single neuron 
should be able to identify the role of each stage in information transfer in terms of the pa- 
rameters relating to the neuron's dendritic structure, its spiking mechanism, etc. Employing 
this reductionist approach, we focus on a important component of neural processing, the 
synapse, and analyze a simple model of a cortical synapse under two different representa- 
tional paradigms. Under the "signal estimation" paradigm, we assume that the input signal 



is linearly encoded in the mean firing rate of a Poisson neuron and the mean-square error 
in the reconstruction of the signal from the post-synaptic voltage quantifies system per- 
formance. From the performance of the optimal linear estimator of the signal, a lower 
bound on the capacity for signal estimation can be computed. Under the "signal detection" 
paradigm, we assume that information is encoded in an all-or-none format and the error in 
deciding whether or not a presynaptic spike occurred by observing the post-synaptic voltage 
quantifies system performance. This is similar to the conventional absent/present(Yes-No) 
decision paradigm used in psychophysics. Performance of the optimal spike detector in 
this case allows us to compute a lower bound on the capacity for signal detection. 
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Figure 1: Schematic block diagram for the signal detection and estimation tasks. The 
synapse is modeled as a binary channel followed by a filter h(t) — at exp(—t/t s ). where 
a is a random variable with probability density, P(a) = a (aa) k ^ 1 exp(—aa) / \k — 1)!. 
The binary channel, (inset, eo = Pr[spontaneous release], t\ = Pr [release failure]) models 
probabilistic vesicle release and h(t) models the variable epsp size observed for cortical 
synapses. n(t) denotes additive post-synaptic voltage noise and is assumed to be Gaussian 
and white over a bandwidth B n . Performance of the optimal linear estimator (Wiener 
Filter) and the optimal spike detector (Matched Filter) quantify synaptic efficacy for signal 
estimation and detection respectively. 



2 The Synaptic Channel 

Synaptic transmission in cortical neurons is known to be highly random though the role of 
this variability in neural computation and coding is still unclear. In central synapses, each 
synaptic bouton contains only a single active release zone, as opposed to the hundreds or 
thousands found at the much more reliable neuromuscular junction. Thus, in response to an 
action potential in the presynaptic terminal at most one vesicle is released (Korn and Faber, 
1991). Moreover, the probability of vesicle release p is known to be generally low (0.1 
to 0.4) from in vitro studies in some vertebrate and invertebrate systems (Stevens, 1994). 
This unreliability is further compounded by the trial-to-trial variability in the amplitude of 
the post-synaptic response to a vesicular release (Bekkers et. al, 1990). In some cases, the 
variance in the size of EPSP is as large as the mean. The empirically measured distribution 
of amplitudes is usually skewed to the right (possibly biased due the inability of measuring 
very small events) and can be modeled by a Gamma distribution. 

In light of the above, we model the synapse as a binary channel cascaded by a random am- 
plitude filter (Fig. 1). The binary channel accounts for the probabilistic vesicle release. e 



and ei denote the probabilities of spontaneous vesicle release and failure respectively. We 
follow the binary channel convention used in digital communications (ei = 1— p), whereas, 
p is more commonly used in neurobiology. The filter h(t) is chosen to correspond to the 
epsp profile of a fast AMPA-like synapse. The amplitude of the filter a is modeled as ran- 
dom variable with density P(a), mean fi a and standard deviation <r a . The CV (standard 
deviation/mean) of the distribution is denoted by CV a . We also assume that additive Gaus- 
sian voltage noise n(t) at the post-synaptic site further corrupts the epsp response. n(t) is 
assumed to white with variance a 2 and a bandwidth B n corresponding to the membrane 
time constant r. One can define an effective signal-to-noise ratio, SNR = E a /N a , given 
by the ratio of the energy in the epsp pulse, Eh = J °° h 2 (t) dt to the noise power spectral 
density, iV = <r 2 /B n . The performance of the synapse depends on the SNR and not on 
the absolute values of Ef t or <x„. In the above model, by regarding synaptic parameters as 
constants, we have tacitly ignored history dependent effects like paired-pulse facilitation, 
vesicle depletion, calcium buffering, etc, which endow the synapse with the nature of a 
sophisticated nonlinear filter (Markram and Tsodyks, 1997). 
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3 Signal Estimation 

Let us assume that the spike train of the presynaptic neuron can be modeled as a doubly 
stochastic Poisson process with a rate A(i) = k(t) * m(t) given as a convolution between 
the stimulus m(t) and a filter k(t). The stimulus is drawn from a probability distribution 
which we assume to be Gaussian. k(t) = exp(—t/r) is a low-pass filter which models 
the phenomenological relationship between a neuron's firing rate and its input current, r 
is chosen to correspond to the membrane time constant. The exact form of k{t) is not 
crucial and the above form is assumed primarily for analytical tractability. The objective is 
to find the optimal estimator of m(t) from the post-synaptic voltage v(t), where optimality 
is in a least-mean square sense. The optimal mean-square estimator is, in general, non- 
linear and reduces to a linear filter only when all the signals and noises are Gaussian. 
However, instead of making this assumption, we restrict ourselves to the analysis of the 
optimal linear estimator, rh(t) = g(t) * v(t), i.e. the filter g(t) which minimizes the 
mean-square error E — ((m(t) — rh(t)) 2 ) where (.) denotes an ensemble average. The 
overall estimation system shown in Fig. 1 can be characterized by an effective continuous 
channel (Fig. 2a) where h{t) = m(t) — m(t) denotes the effective reconstruction noise. 
System performance can be quantified by E, the lower E, the better the synapse at signal 
transmission. The expression for the optimal filter (Wiener filter) in the frequency domain is 
g(u>) = S mv (— u})/S vv (ijj) where S mv (u)) is the cross-spectral density (Fourier transform 
of the cross-correlation R mv ) of m(t) and s(t) and S vv (co) is the power spectral density of 
v(t). The minimum mean-square error is given by , E = v^— J s | S mv (u>) | 2 / S vv (uj) cku. 
The set S = {uo \ S vv (u) ^ 0} is called the support of S vv (w). 



Another measure of system performance is the mutual information rate 7(m; v) between 
m(t) and v(t), defined as the rate of information transmitted by v(t) about s(t). By the 
Data Processing inequality (Cover 1991), I(m, v) > I(m, rh). A lower bound of I(m, rh) 

and thus of 7(m; v) is given by the simple expression = \ j s \og 2 [ S g! m ^ ] dm (units 
of bits/sec). The lower bound is achieved when n(t) is Gaussian and is independent of 
m(t). Since the spike train s(i) = $[t ~ h) is a Poisson_process with rate k(t) * m(t), 
its power spectrum is given by the expression, S ss (oj) = A+ | K(cv) | 2 S mm (uj) where 
A is the mean firing rate. We assume that the mean (/x m ) and variance (cr 2 ^) of m(t) are 
chosen such that the probability that A(i) < is negligible^ The vesicle release process 
is the spike train gated by the binary channel and so it is also a Poisson process with rate 
(1 — ei)A(t). Since v(t) = a,ih(t — ti) + n(t) is a filtered Poisson process, its power 
spectral density is given by S vv (uj) =| H(w) | 2 {(^+o^)(l-ei)A+/x 2 ! (l-e 1 ) 2 | K(oj) | 2 
S mm (u))} + S nn (oj). The cross-spectral density is given by the expression S vm (u>) — 
(1 — €i)/j,a,S m m(u)H(oj)K(oj). This allows us to write the mean-square error as, 



Js -W/M + S mm {u) + S e ff(u) 

\{i + cv2) s nn (u) 



(1 - ei) | K{u) p ' ' (1 - ei ) 2 M 2 | H(w) | 2 | K(w) | 2 

Thus, the power spectral density of n(t) is given by Snh = ^e//(^) + S e ff(ui). Notice 
that if K(u>) — > oo, E — > /.e. perfect reconstruction takes place in the limit of high 
firing rates. For the parameter values chosen, S e ff(u) -C A e //(w), and can be ignored. 
Consequently, signal estimation is shot noise limited and synaptic variability increases shot 
noise by a factor N syn = (1 + CV^)/(l - ei). For CV a = 0.6 and e x = 0.6, N syn = 3.4, 
and for CV a = 1 and ei = 0.6, iV sa „ = 5. If m(t) is chosen to be white, band-limited to 
B m Hz, closed-form expressions for E and In, can be obtained. The expression for In, is 
tedious and provides little insight and so we present only the expression for E below. 

Eft,!*) = o* [1 - l-ta n -\-^)] 



E is a monotonic function of 7 (decreasing) and By (increasing). 7 can be considered as 
the effective number of spikes available per unit signal bandwidth and Bt is the ratio of 
the signal bandwidth and the neuron bandwidth. Plots of normalized reconstruction error 
E r = EJ £r 2 ra and In, versus mean firing rate (A) for different values of signal bandwidth B m 
are shown in Fig. 3a and Fig. 3b respectively. Observe that In, (bits/sec) is insensitive to B m 
for firing rates upto 200Hz because the decrease in quality of estimation (E increases with 
B m ) is compensated by an increase in the number of independent samples (2B m ) available 
per second. This phenomenon is characteristic of systems operating in the low SNR regime. 

has the generic form, 7^ = B log(l + S/(NB)), where B, S and TV denote signal 
bandwidth, signal power and noise power respectively. For low SNR, I w B S/ (NB) = 
S/N, is independent of B. So one can argue that, for our choice of parameters, a single 
synapse is a low SNR system. The analysis generalizes very easily to the case of multiple 
synapses where all are driven by the same signal s(t). (Manwani and Koch, in preparation). 
However, instead of presenting the rigorous analysis, we appeal to the intuition gained from 
the single synapse case. Since a single synapse can be regarded as a shot noise source, 
n parallel synapses can be treated as n parallel noise sources. Let us make the plausible 

1 We choose fj, m and a m so that A — 3a\ (std of A) so that Prob[A(t) < 0] < 0.01. 



assumption that these noises are uncorrelated. If optimal estimation is carried out separately 
for each synapse and the estimates are combined optimally, the effective noise variance 
is given by the harmonic mean of the individual variances i.e. l/a^ e fj — J2i l/ a ni- 
However, if the noises are added first and optimal estimation is carried out with respect 
to the sum, the effective noise variance is given by the arithmetic mean of the individual 
J2i a ni/ n2 - ^ we assume that all synapses are similar so that 



variances, i.e. a. 



neff 

a 'ni — °' 2 ' a neff = ® l l n - Plots of E r and for the case of 5 identical synapses are 
shown in Fig. 3c and Fig. 3d respectively. Notice that increases with B m suggesting 
that the system is no longer in the low SNR regime. Thus, though a single synapse has very 
low capacity, a small amount of redundancy causes a considerable increase in performance. 
This is consistent with the fact the in the low SNR regime, I increases linearly with SNR, 
consequently, linearly with n, the number of synapses. 
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Figure 3: E r and lib vs. mean firing rate (A) for n = 1 [(a) and (b)] and n = 5 [(c) and (d)] identical 
synapses respectively (different values of B m ) for signal estimation. Parameter values are ei = 0.6, 
e = 0, CV a = 0.6, t 3 = 0.5 msec, r = 10msec, a„ = 0.1 mV, B n = 100 Hz. 



4 Signal Detection 

The goal in signal detection is to decide which member from a finite set of signals was 
generated by a source, on the basis of measurements related to the output only in a statistical 
sense. Our example corresponds to its simplest case, that of binary detection. The objective 
is to derive an optimal spike detector based on the post-synaptic voltage in a given time 
interval. The criterion of optimality is minimum probability of error (P e ). A false alarm 



(FA) error occurs when a spike is falsely detected even when no presynaptic spike occurs 
and a miss error (M) occurs when a spike fails to be detected. The probabilities of the errors 
are denoted by Pj and P m respectively. Thus, P e = (1 — p Q ) Pj +p Q P m where p Q denotes 
the a priori probability of a spike occurrence. Let X and Y be binary variables denoting 
spike occurrence and the decision respectively. Thus, X — 1 if a spike occurred else X — 
0. Similarly, Y = 1 expresses the decision that a spike occurred. The posterior likelihood 
ratio is defined as C(v) = Pr(v \ X = 1)/Pr(v \ X = 0) and the prior likelihood as 
C = (1 — p )/p . The optimal spike detector employs the well-known likelihood ratio 
test,"If£(w) > C Y = lelseY = 0". WhenX = l,v(t) =ah(t)+n(t) elseu(t) = n(t). 
Since a is a random variable, C(v) = (J Pr(v \ X = l;a) P(a) da)/Pr(v \ X = 0). If 
the noise n{t) is Gaussian and white, it can be shown that the optimal decision rule reduces 
to a matched filter^ i.e. if the correlation, r between v(t) and h(t) exceeds a particular 
threshold (denoted by rj), Y = 1 else Y = 0. The overall decision system shown in 
Fig. 1 can be treated as effective binary channel (Fig. 2b). The system performance can 
be quantified either by P e or I(X; Y), the mutual information between the binary random 
variables, X and Y. Note that even when n(t) — (SNR — oo), P e ^ due to the 
unreliability of vesicular release. Let P* denote the probability of error when SNR = oo. 
If eo = 0, P* = p ei is the minimum possible detection error. Let PJ and P° n denote FA 
and M errors when the release is ideal (ei = 0, eo = 0)- It can be shown that 

P e = P: + P° m [p {l - ex) - (1 -p o )e ] + P° f [{l-p )(l - e )- Po e 1 ] 
p f = p; , P m = P 7 ° n + - P° + PJ) 

Both P^ and P^ depend on r\. The optimal value of r\ is chosen such that P e is minimized. 
In general, PJ and P^ can not be expressed in closed-form and the optimal r\ is found using 
the graphical ROC analysis procedure. If we normalize a such that /i a — 1, Pj? and P^ can 
be parametrically expressed in terms of a normalized threshold ry*, Pj? = 0.5 [1 — Erf(rj*)}, 
P° m = 0.5[1 + J °° Erf(rf - VSNR a) P(a) da}. I{X; Y) can be computed using the 
formula for the mutual information for a binary channel, / = 7i(p (1 — P m ) + (1 — 
Po) Pf)-PoU{P m ) - (l-p )W(P/) where U{x) = -xlog 2 (x) - (l-;r)log 2 (l-x) is 
the binary entropy function. The analysis can be generalized to the case of n synapses but 
the expressions involve ri-dimensional integrals which need to be evaluated numerically. 
The Central Limit Theorem can be used to simplify the case of very large n. Plots of 
P e and I(X; Y) versus n for different values of SNR (l,10,oo) for the case of identical 
synapses are shown in Fig. 4a and Fig. 4b respectively. Yet again, we observe the poor 
performance of a single synapse and the substantial improvement due to redundancy. The 
linear increase of / with n is similar to the result obtained for signal estimation. 

5 Conclusions 

We find that a single synapse is rather ineffective as a communication device but with 
a little redundancy neuronal communication can be made much more robust. Infact, a 
single synapse can be considered as a low SNR device, while 5 independent synapses 
in parallel approach a high SNR system. This is consistently echoed in the results for 
signal estimation and signal detection. The values of information rates we obtain are very 
small compared to numbers obtained from some peripheral sensory neurons (Rieke et. al, 
1996). This could be due to an over-conservative choice of parameter values on our part 
or could argue for the preponderance of redundancy in neural systems. What we have 
presented above are preliminary results of work in progress and so the path ahead is much 

2 For deterministic a, the result is well-known, but even if a is a one-sided random variable, the 
matched filter can be shown to be optimal. 
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P e (a) and (b) vs. the number of synapses, n, (different values of SNR) for signal detection. 
SNR = Inf. corresponds to no post-synaptic voltage noise. All the synapses are assumed to be 
identical. Parameter values are p a = 0.5, ei = 0.6, eo = 0, CV a = 0.6, t s = 0.5 msec, r = 10 msec, 
a n = 0.1 mV, B n = 100 Hz. 



longer than the distance we have covered so far. To the best of our knowledge, analysis 
of distinct individual components of a neuron from an communications standpoint has not 
been carried out before. 
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